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Summary : 

Even  for  Reynolds  numbers  not  small,  the  dynamic 
effects  including  the  dynamic  velocity  terms  and  the 
pressure  terms  enumerated  in  the  Navier  Stokes  equation 
exert  small  average  effects  on  the  turbulence  features 
when  compared  with  the  effects  of  viscosity. 

This  new  finding  points  to  the  linearization 
of  the  turbulence  problem.  Non-linear  effects  of  the 
statistical  deviations  from  the  random  distribution, 
however,  are  not  insignif icent . 

The  theory  is  applied  to  the  determination  of  the  ratio 
of  the  two  intensities  of  turbulence  mixing,  the  one  being  the 
shear  mixing  and  the  other  the  energy  diffusion  mixing. 
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I.  Introduction, 

The  turbulent  fluid  motion  presents  a  problem  in  statistics. 

It  must  be  the  aim  of  any  real  turbulence  theory  to  sift  experi¬ 
mental  and  theoretical  knowledge  to  the  utmost,  so  that  we  may  be 
able  to  predict  a  small  number  of  relevant  averages  by  the  use  of 
as  small  a  number  of  averaged  variables  as  possible. 

We  have  shown  in  Ref .  ( m)  that  even  the  restriction  to  three 
such  variables  affords  a  good  insight  into  the  general  features  of 
turbulent  fluid  motion.  This  encourages  us  to  go  forward  in  the 
chosen  direction.  The  present  paper  deals  with  refinements  of  the 
mathematical  procedure.  It  brings  our  approach  in  line  with  the 
branch  of  the  current  literature  generally  summed  up  as  "statisti¬ 
cal  theory  of  turbulence."  This  is  not  difficult  in  principle  be¬ 
cause  an  average  lump  flow  is  not  quite  unrelated  to  the  components 
of  the  velocity  correlation  dyadic.  The  lump  picture  and  the  abstract 
correlation  view  are  in  harmony  with  each  other.  The  correlations 
appeal  more  to  the  abstract  mathematician.  The  lumps  seem  to  offer 
more  hope  of  progress  to  the  physicist  who  is  aware  of  the  danger 
of  pitfalls  of  pure  abstract  reasohing. 

Time  and  space  forbids  -us  to  present  a  self-contained  account 
of  the  current  statistical  approach.  The  following  is  accordingly 
directed  to  readers  who  are  reasonably  acquainted  with  at  least  the 
two  pionefer  papers  on  the  subject,  by  G.I.  Taylor  and  Theodore  von 
Karmah.  Ref.  (t)  and  ( k) .  A  relatively  recent  account  by  Agostini 
may  also  be  -consulted  with  profit.  It  gives  a  good  idea  about  the 
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present  status  .of  the  theory,  about  its  satisfactory  aspects  as  well 
as  its  weak  ones,  'An  excellent  translation  of  Agostini's  paper  has 
been  made  by  J.  Vanier,  and  is  readily  available.  Ref.  (a), 

-T-I .  The  postulates  of  the  statistical-’  theory . 

All  postulates  employed  in  the  statistical  theory  of  turbulence 
can  be  classified  under  the  heading  of  three  different  principles. 

They  flow  either  from  the  physical  theory,  or  from  an  attempt  to  sepa¬ 
rate  a  turbulence  motion  from  a  basic  motion,  or  lastly  from  some  kind 
of  artificial  specialization. 

The  physical  theory  by  itself  is  not  uncertain  nor  controver¬ 
sial.  It  is  exhaustively  laid  down  in  the  Navier  Stokes  equations 
in  combination  with  the  equation  of  continuity.  It  Is  Implied  that 
the  solutions  deoend  also  and  in  a  mosit  pronounced  way  on  the  boundary 
conditions  clearly  necessary  or  pierely  assumed.  The  Navier  Stokes 
equations  are  ordinarily  stated  with  reference  to  incompressible  fluids 
only.  When  adopted  In  that  special  form,  it  Is  out  of  place  to  specu¬ 
late  again  whether  the  fluid  Is  assumed  to  be  Incompressible  or  not. 

smooth 

The  division  of  the  actual  flow  Into  a  basics  flow  and  a  turbu¬ 
lence  flow  is  quite  generally  resorted  to.  But  it  is  not  always 
realized  that  this  division  ip  not  unique.  It  will  not  do  just  to 
say  that  the  basic  flow  is  the  average,  and  the  turbulence  flow  is  what 
remains  over  and  above  the  average.  That  is  not  precise  enough.  The 
division  has  also  to  Include  the  boundary  conditions.  Also  the 
averaging  process  can  be  carried  out  in  more  than  one  manner  and  to 
different  degree  and  intensity.  There  is  in  fact  no  sharp  demarca- 

9 

tlo’n  line  between  the  turbulence  flow  and  the  basic  flow.  If  we 
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c&rry  the  averaging  too  far  all  motion  becomes  turbulence  motion.  If 
hot  far  enough,  the  basic  flow  remains  too  involved. 

Specialization,  such  as  to  homogeneous  flows  and  to  Isotropic 
turbulence,  is  resorted  to  for  dividing  the  difficulties  of  the 
problem.  It  is  attempted  to  investigate  in  that  fashion  different 
aspects  of  turbulence  by  themselves.  This  is  done  with  the  expecta¬ 
tion  that  later  these  different  effects  may  be  combined  or  superposed 
to  each  other  in  some  simple  way.  This  subdivision  and  the  postulates 
necessary  to  effect  it  are  again  not  unidjue  nor  precise,  but  they  are 
particularly  tricky  and  readily  lead  to  gapi  in  the  logical  procedure. 

An  orderly  progress  in  the  presentlly  discussed  theory  requires 
that  the  postulates  necessary  for  applying  these  three  principles 
be  Introduced  gradually  and  with  greatest  care.  The  decisions  should 
be  delayed  until  they  become  immediately  necessary  In  each  individual 
case..  That  particular  case  serves  then  for  illustrating  and  explain¬ 
ing  the  scope  of  the  postulate.  It  leads  only  to  confusion  if  spme 
so-called  definition  is  accepted  as  a  dogma  to  serve  as  the  basis  of 
the  consequent  research.  The  mathematical  procedure  includes  then 
"proving"  that  certain  relations  hold  by  trying  to  demonstrate  that 
they  conform  with  the  perhaps  ambiguously  stated  or  prematurely 
adopted  definitions. 

It  may  seem  strange  that  such  reserved  attitude  should  also  be 
necessary  with  regard  to  the  physical  theory,  which  stands  establish¬ 
ed  certain  and  rigid.  But  the  boundary  conditions  are  a  matter  of 
discretion,  and  cannot  be  established  without  regard  to  physics.  Navier 
Stokes  Stokes  equations  do  not  determine  the  boundary  conditions. 
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Said  equations  are  precise  when  applied  to  specific  solutions.  But 
there  is  no  solution  to  apply  them  to.  It  is  the  very  nature  of  the 
turbulence  problem  that  the  precise  solution  must  forever  remain  un¬ 
known  and  unstated.  The  application  of  the  equations  in  such  cases 
requires  mental  steps  not  apparent  from  the  equation's.  The  formu¬ 
lated  equations  stand  for  broad  principles,  not  for  a  narrow  numer¬ 
ical  procedure.  A  full  grasp  of  all  implications  of  the  equations 
has  to  be  constantly  resorted  to. 

All  applied  statistical  theory  is  attended  with  some  lack  of 
rigour.  Thus  for  instance,  the  standard  or  Gaussian  distribution  is 
often  used  even  if  it  16  not  the  true  distribution.  The  theory  deals 
with  incomplete,  scant  and  uncertain  data,  and  conclusions  regarding 
a  shape  have,  so  to  dpeak,  to  be  drawn  from  its  vague  shadow.  Under 
such  circumstances  we  cannot  afford  over-fine  mathematical  distinct¬ 
ions  to  stand  in  the  way  ot  using  quite  workable  and  practical  propo¬ 
sitions,  improvable  as  they  may  be.  The  door  to  ingenuous  and  shrewd 
guesses  must  not  be  closed.  Simplifying  assumptions  should  not  be 
rejected  for  lack  of  mathematical  rigours  unless  it  is  demonstrated 
that  they  lead  to  unacceptable  conclusions.  In  such  manner  was  the 
kinetic  theory  of  gases  developed.  The  turbulence  problem  is  a  much 
humbler  problem,  having  no  significance  by  Itself  but  only  having  a 
significance  on  account  of  engineering  needs.  Certainly,  what  is  good 
enough  for  fundamental  physics  is  good  enough  for  turbulence  research. 


Ill 
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IXI.  The  correlation  flow. 

The  correlation  flow  is  obtained  by  a  scalar'  multiplication 
of  the  correlation  dyadic  (tensor  of  second  rank) 

V  =  'O  ■  U'i  V  /to'V 

R  is  the  correlation  dyadic  u;  v,  where  the  semicolon  has  no  signi¬ 
ficance  beyond  that  of  the  empty  space  ordiparily  employed  for 
indicating  tensor  multiplication.  The  symbol  "1"  denotes  a  unit  vector. 
An  entire  tensor  field  R  is  associated  with  each  point  of  the  space. 
Each  field  is  in  the  associated  space,  and  only  v  but  not  u  varies 
in  this  space.  All  variables  depend  besides  on  time, all  averages 
are  time  averages.  The  divergence  is  zero 

^7‘V  -  O 

because  averaging  leaves  absence  of  divergence  unchanged. 

It  is  always  postulated  that  the  turbulence  is  reasonably 
homogeneous  so  that  the  changes  of  the  averages  through  the  territory 
of  one  lump  are  negligibly  small  when  compared  with  the  averages 
themselves.  This  postulate  is  often  unrealistic. 

The  present  discussion  is  restricted  throughout  to  homogeneous 
and  incompressible  fluids.  The  correlation  flow  is  the  average  of 
the  turbulence  motion  at  one  point  provided  the  average  is  only  taken 
while  the  velocity  at  that  point  is  reasonably  near  to  one  predeter¬ 
mined  direction.  That  Is  certainly  so  for  isotropic  turbulence.  We 
see  or  assume  that  the  correlation  flow  possesses  a  number  of  symmetric 
features.  It  is  ordinarily  assumed  to  possess  an  axis  of  symmetry 
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parallel  to  the  central  velocity.  This  does  not  mean  yet  isotropic 
turbulence.  We  shall  see  or  assume  that  the  correlation  flow  possess¬ 
es  a  number  of  symmetric  features.  It  is  ordinarily  assumed  to  possess 
an  axis  of  symmetry  parallel  to  the  central  velocity.  This  does  not 
mean  yet  isotropic  turbulence.  Only  if  the  correlation  flow  is  the 
same  for  all  directions  is  the  turbulence  isotropic.  We  hope  to  be  able 
to  reach  tangible  results  even  in  the  case  that  the  turbulence  is  only 
approximately  ad.  merely  to  some  degree  Isotropic,  by  using  suitable 
averages  of  the  correlation  flows  in  different  directions. 

The  correlation  flows  reach  theorically  up  to  the  infinity. 

They  converge  rapidly  to  zero,  as  we  shall  discuss  in  fuller  detail 
in  Section  5-  Still  it  is  impossible  to  assign  a  sharp  boundary  up  to 
which  the  flows  actively  reach.  A  similar  situation  was  already  faced 
with  our  Hill  vortices  in  Ref .  (m),  and  is  quite  common  in  applied 
physics . 

There  is  a  separate  field  or  space  function  of  correlation 
dyadics  and  therefore  also  one  of  correlation  flows  for  each  point. 

In  that  respect  the  correlation  flow  is  different  from  the  lump  flow. 
Also  lumps  have  a  physical  reality,  but  correlations  are  merely  prob¬ 
ability  weight  functions.  They  change  with  the  change  of  the  entire 
flow  but  they  have  no  individual  history,  because  they  have  no  in¬ 
dividual  physical  existence.  Thus,  in  a  quasi-steady  turbulent  flow 
the  correlations  and  therefore  the  correlation  flows  do  not  change 

The  individual  lumps  change  each  by  themselves  even  their  average 
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never  changes.  Th'e  lump  flow  has  always  a  history  of  change.  It  is 
the  plan  of  our  liiquiry  to  study  the  average  history  of  the  individual 
lump  flows  as  already  done  in  Ref.  (m). 

The  average  lump  flows  are  closely  related  to  the  core  or  inner 
portion  of  the  correlation  flow.  The  two  are  not  strictly  identical 
but  they  suggest  each  other.  The  lump  flows,  even  the  average  lump 
flows,  are  a  physical  reality  and  it  is  not  obvious  that  a  certain 
correlation  flow  corresponds  to  one  lump  flow  only.  The  study  of 
such  detail  questions  is  however  contrary  to  the  plan  of  a  statistical 
approach.  We  shall  never  know  either  of  the  two  flows  with  any  pre¬ 
cision.  The  distinction  does  not  appear  relevant  at  this  state  of  the 
inquiry . 

We  plan  to  disregard  the  distinctions  between  the  core  of 
the  correlation  flow  and'  the  lump  flow  ,  except  for  a  proportionality 
factor  for  the  velocity  magnitudes.  The  strength  of  the  lump  flow 
in  relation  to  the  strength  or  energy  density  of  the  turbulence  de¬ 
pends  on  the  average  lump  distance.  We  plan  to  present  an  inquiry 
about  that  in  a  later  paper.  The  present  paper  deals  with  the 
supremacy  of  the  viscosity  only. 
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IV.  Correlation  dynamics 

In  current  literature  on  the  statistical  theory  of  turbulence 
the  term  "dynamics"  designates  the  study  of  the  time  changes  of  the 
correlation  dyadic  of  pairs  of  simultaneous  velocities.  If  the  flow 
is  quasi-steady,  there  is  not  any  such  change.  But  the  theory  is 
mostly  applied  to  isotropic  turbulence.  Then  the  turbulence  always 
decays  and  the  correlation  dyadic  changes. 

The  fluid  has  to  rely  on  pressure  effects  for  the  preservation 
of  density.  The  current  literature  notes  that  nevertheless,  for 
isotropic  turbulence,  the  pressure  terms  of  the  Navler  Stokes  equation 
cancel  out  by  averaging,  and  leave  no  trace  in  the  equation  describing 
the  history  of  the  correlation.  The  correlation  dyadic  undergoes 
changes  Just  as  if  there  were  no  pressures  acting.  This  propo¬ 
sition  is  ordinarily  accepted  on  the  authority  of  von  Karman, 

Bef.(k).  The  proof  or  demonstration  given  there  seems  to  be  not 
generally  accepted.  Most  authors  do  not  discuss  this  point.  Agostini, 
Hef.  (a),  page  69,  discusses  the  disappearance  of  the  scalar  product 
of  the  velocity  and  the  pressure  gradient,  not  their  dyadic  product. 
Both  products  show  equal  response  or  non-response  to  the  pressure. 
Agostini  does  not  accept  the  proposition.  A  few  pages  later,  Hef. (a) 
p.  75,  he  says:  "When  assuming  the  disappearance  of  the  pressure  term 
to  have  proved."  In  effect  he  repudiates  with  these  words  the  proof 
he  has  himself  offered  and  rejects  the  different  proof  of  ' von  Karman 
which  he  must  have  studied  carefully.  The  present  author,  while 
realizing  that  Hef.  ( k)  Section  7  is  perhaps  not  quite  as  lucid  as  it 
could  be  in  view  of  the  simplicity  and  directness  of  the  reasoning 
presented,  still  does  not  find  any  error  or  fault  either  in  the 
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postulates  employed  nor  In  the  rather  primitive  and  elementary  steps 
taken.  The  pressure  tterms  do  disappear.  It  is  however  noted  that  the 
demonstration!  in  said  Section  7  holds  equally  for  any  other  scalar 
variable  subject  to  the  condition  of  statistical  Isotropy.  Nothing 
is  introduced  in  Section  7  having  reference  to  the  scalar  rWr' 
representing  the  pressure.  This  widens  the  scope  of  application  of 
the  result  of  Section  7  in  that  it  applies  to  all  Isotropic  scalar 
variables,  but  \t  also  narrows  it  unnecessarily  in  that  it  restricts 
the  disappearance  of  the  correlation  dyadic  to  a  product  of  the 
velocity  and  of  a  vector  belonging  to  an  irrotational  Victor  field. 

f 

The  final  outcome  appears  in  quite  a  different  light  when 
relying  more  on  physical  grounds  and  less  on  mere  geometry  and 
postulates.  The  pressure  terms  are  then  seen  to  average  out  without 
having  to  rely  on  pressure  gradients  being  derived  from  a  scalar 
potential.  They  average  out  for  a  much  broader  reason,  namely  because 
the  pressure  is  an  even  function  of  the  velocity.  If  all  turbulence 
velocities  are  permitted  to  be  reversed,  the  pressure  distribution 
would  remain  what  it  was  before,  without  change  of  sign.  It  is  there¬ 
fore  reasonable  to  assume  that  for  each  pressure  value  all  velocities 
occur  as  often  in  one  direction  and  magnitude  as  in  the  opposite 
direction  and  the  same  magnitude.  This  leads  immediately  to  the 
disappearance  of  the  pressure'-velocity  correlation  gradient.  Neither 
a  potential  of  the  pressure  force,  nob  Isotropy  of  the  turbulence  motion 
is  necessary  for  proof.  Now,  this  apparently  new  mode  of  reasoning 
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is  of  larger  significance  than  just  furnishing  a  new  proof  for  an 
old  and 'long  accepted  proposition.  It  widens  the  scope  of  appli¬ 
cation  to  less  restricted  turbulence  and  to  other  terms  of  the 
equation.  It  directs  our  attention  to  the  remaining  dynamic  pres¬ 
sure  terms  of  the  Navier  Stokes  equation.  They  can  be  written 

V-  (  v;  v)  and  they  are  also  an  even  function  of  the  velocity. 
Exactly  the  same  reasoning  that  holds  for  the  pressure  terms  applies 
also  to  the  dynamic  velocity  terms.  These  terms  will  also  cancel 
out  on  the  average  for  the  history  of  the  correlation  dyadic  and  of 
correlation  flow.  That  puts  the  theory  in  an  entirely  new  light. 
Only  the  viscosity  terms  remain,  governing'  the  correlation  dyadic 

^  and  the  correlation  flow,  even  if  the 
Reynolds  number  is  not  particularly  small. 

Retaining  merely  the  viscosity  term  was  already  considered  by 
von  Karman,  Ref.  (k)  Section  10,  but  expressly  restricted  to  "small 
Reynolds  numbers."  We  see  now  that  the  supremacy  of  the  viscosity 
is  quite  general.  Agostini,  Ref.  (a)  Section  20,  also  discusses 
this  discounting  of  the  dynamic  velocity  terms,  but  with  apology 
and  little  faith.  He  defends  it  by  stating  that  it  is  mathematic¬ 
ally  convenient  but  physically  quite  difficult  to  justify.  He  over¬ 
looked  the  physical  justification  presented  in  the  present  section 
believing  rather  that  there,  was  none. 


The  correlation  distribution  is  thus  seen  to  be  subject  to 
a  linear  differential  equation.  This  situation  must  however  not 
be  accepted  without  reservation.  What  happens  in  the  turbulence 
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flow  does  not  depend  merely  on  the  probability  distribution  but  also 
on  the  deviations  from  that  probability.  Secondary  deviation  effects 
must  not  be  overlooked.  Section  7  of  this  paper  illustrates  this 
remark. 

V.  The  correlation  flow  at  large  distances. 

It  is  necessary  to  clarify  the  meaning  of  the  total  momentum 
integral,  that  is  the ' space  integral  of  the  correlation  velocity 
vector, 

5-1 

the  integration  extending  through  the  entire  space.  Also  that  of 
the  space  integral  for  the  total  polar  moment  of  inertia  of  the 
correlation  velocity, 

5-2  V  cd  S' 

where  v  denotes  the  velocity,  r  denotes  the  radial  distance  from 
the  center,  and  dS  denotes  the  space  elements. 

With  one  Hill  vortex  or  with  a  similar  correlation  flow  the 
velocities  decline  as  the  cube  of  the  radial  distance  r.  Hence  the 
momentum  integral  (5-1)  does  not  converge.  Its  value  depends  on 
the  order  in  which  the  integration  is  carried  out,  and  on  the 
specification  for  the  infinitely  far  boundary.  Integral  (5-2) 
would  even  less  converge. 

We  have  to  look  deeper  into  this  question  of  convergence.  We 
must  call  to  mind  that  we  are  facing  not  one  single  lump  but  an 
infinite  assembly  of  lumps,  moving  in  some  random  fashion  about  in 
the  fluid.  Their  distribution  is  not  entirely  random  in  that  there  is 
a  bias  or  correlation  between  the  motion  of  adjacent  lumps.  This  is 
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the  case  by  reason  of  the  physical  laws  governing  what  happened  in 
the  past.  Individual  lumps  have  only  a  short  life  span,  and  they 
are  all  created  by  the  action  of  the  same  pressure  distribution. 

The  momentum  of  the  fluid  is  not  changed  but  merely  redistributed. 

In  consequence  each  lump  traveling  in  one  direction  is  surrounded 
by  a  group  of  lumps  having  on  the  average  a  total  momentum  in  the 
opposite  direction.  Hence,  the  correlation  spreads  out  far  beyond 
the  territory  of  the  lump.  Also  this  momentum  is  evenly  distributed 
around  the  first  mentioned  lump.  There  is  no  upper  limit  for  the 
statistical  symmetry  thus  established  *  In  consequence  ,  the  corre¬ 
lation  velocity  declines  much  faster  than  inverse  to  the  third  power 

of  the  distance.  It  would  at  least  decline  with  the  fifth  power,- 

which 

but  there  is  no  upcer  limit  of  the  exponent  m  byA the  velocity  may 

-Kti 

decline  with  the  distance.  An  exponential  law  g  seems  appro¬ 
priate  and  acceptable . 

This  must  not  be  taken  too  literally.  The  sense  of  reality 
must  be  preserved.  The  higher  the  power,  the  slower  the  convergence, 
and  the  larger  the  deviation  from  the  mean.  It  will  require  a  longer 
and  longer  time  to  obtain  the  average  and  larger  space  too.  The 
average  space  integrals  become  very  unreal  as  the  power  m  goes  up. 

It  is  Impossible  to  determine  them  experimentally.  But  for  the 
following  we  do  not  have  to  go  very  high  with  the  exponent  m.  In 
consequence  of  this  physical  reasoning  we  feel  justified  in  postulating 
thgt  all  correlation  integral#  occurring  converge  at  large  distance 
from  the  origin. 
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VI.  On  the  preservation  of  the  lump  flow. 

Additional  comment  regarding  the  minor  part  played  by  the 

inertia  effects  in  the  turbulence'  theory  seems  to  be  called  for. 

vIt  J,s  clear  that  each  individual'  lump  is  attended  by  a  pressure 

distribution  apd  by  inertia  effects.  We  have  already  pointed  out 

in  Ejection  4  why  these  inertia  effects  are  secondary  with  regard 

"*  ■  a 

to  average  effects.  ■  We  would  like  to  show  this  ini  more  tangible 
way  using  the  physical  picture  of  one  lump  undergoing  changes. 

If  a  lump,  under  the  action  of  inertia  only,  preserves  its 
compactness,  the  flow  picture  will,  to  some  extent,  be  more  or  less 
periodical,  and  therefore  its  time  average  will  be  constant.  The 
basic  question  is  accordingly  whether  a  compact  lump  flow,  of  the 
type  considered,  which  is  a  generalized  vortex  ring  in  a  non-viscous 
fluid  similar  to  the  Hill  vortex  of  Ref.  ( m)  remains  compact.  The 
opposite  would  be  that  the  vortices  spread  out  to  become  distributed 
through  a  larger  and  larger  territory,  i.e.  in  the  absence  of  viscosity. 
yUhe  vortex  rings  of  the  correlation  flow  occur  equally  in  both 
directions  but  the  vortex  rings  of  the  lumps  are  predominantly  in  one 
direction.  The  vorticlty  distribution  through  the  outer  territory 
of  the  correlation  flow  can  not  count  for  the  Itftnp-,.  because  that 
consist#  of  many  lumps  rather  than  of  an  even  and  gradual  distribution 
of  vorticity.  Now  in  absence  of  viscosity  effects,  the  total  energy 
of  the  lump  flow  must  remain  constant.  The  total  kinetic  energy  of 
the  flow  depends  on  the  compactness  of  the  vortex  distribution. 

The  wider  the  same  vortex  lines  are  scattered  and  distributed,  the 
smaller  is  the  kinetic  energy  of  the  flow.  Ref  (1)  Section  153 ■ 

It  follows  therefore  from  the 
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energy  principle  together'  with  the  fact  that  predominately  only 

f 

vortex  rings  of  the  same  direction  are  involved,  that  the  lump  flow 
must  remain  compact,  -'-t  is  therefore  statistically  self-preserving;, 
except  for  viscosity  effects. 

VII.  Diffusion 

We  are  interested  in  the  diffusion  of  a  quantity  that  is  con¬ 
centrated  in  the  vicinity  of  a  point  so  that  it  converges  towards  zero 
within  reasonable  distance  fro*  that  point,  as  described  in  Section  5* 
We  have  to  distinguish  between  the  case  that  the  space  integral  of 
the  quantity  is  zero  ?-l  s  -  o 


and  when  it  IS  not.  If  the  space  integral  Is  zero,  the  moment  of  inertia 


of  the  quantity  is  constant 


o/S 


7-2 J/l>  ^  ^  ~  const 

This  follows  from  applying  the  diffusion  equation 

”  £  r  -  y  v-vr 

to(7-2.)and  applying  Green's  theorem,  two  times  in  succession. 

The  surface  integrals  vanish  in  consequence  of  the  assumptions  made 
regarding  the  decline  of  the  quantity  with  increasing  distance, 
Section  5-  Hence  we  obtain  (7-2).  More  general,  the  first  moment 
computed  with  the  (2n)th  power  of  r  not  being  zero  will  remain 
constant.  It  is  practically  sufficient  for  obtaining  (7-2),  that 


(7-1)  be  very  small,  that  is  that  the  Reynolds  number 


(^/V/y 


be  large  where  L  and  t  denote  respectively  the  largest  distance  and 
the  largest  time  span  Involved  in  the  experiment. 
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We  have  to  difetlnguish  between  the  moment  o.f  Inertia  of  the 
quantity  being  finite  or  ?ero .  There  exist  type  preserving  distributions 
of  finite  magnitude  of  the  monjent  of  inertia.  They  are  the  fundamental 
or  dominating  term  of  a  combination  of  type  perservlng  functions. 

They  are  tjhe  ones  having  a  finite  moment  of  Inertia.  The  other  ones 
have  the  moment  of  inertia  zero. 

These  distributions,  type -perservlng  under  the  action  of  viscosity 
only,  have  already  been  introduced  into  the  theory  in  Ref.  ( k) . 

However,  they  were  only  Introduced  fob  small  Reynolds  numbers . 

We  see  now  that  their  field  of  usefulness  is  much  broader.  This 
statistically  determined  effect  of  viscosity  brings  this  about.  For 
the  time  being,  and  quite  possibly  for  all  time,  it  is  necessary 
to  restrict  ourselves  to  using  one  prototype  flow  to  serve  Once  for 
all  as  correlation  flow,  and  also  by  way  of  its  inner  part,  as  pro¬ 
totype  lump  flow.  Going  into  the  details  of  more  than  one  lump 
flow  seems  to  be  quite  at  variance  with  the  policy  of  a  statistical 
approach.  The  orototype  flow  would  take  the  place  of  the  Hill  vortex 
in  Ref.  (rr)  . 

The  dominant  influence  of  the  viscosity  demonstrated  in  Section 
4  calls  distinctly  for  the  choice  of  a  prototype  flow  self-preserving 
in  type  under  viscosity  effects,  not  under  dynamic  effect  as  the' 

Hill  vortex  is.  The  first  choice  would  naturally  be  the  above  funda¬ 
mental  self-preserving  flow-,  having  a  finite  moment  of  the  velocity. 

This  flow  has  also  already  appealed  to  the  authors  of  Ref.  ( k)  and 
(a)  and  has  there, at  least  been  alluded  to.  This  tends  to  show  that 
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the  supremacy  of  the  turbulence  was  intuitively  felt  for  many  years. 

But  it  was  not  distinctly 'reoogpized  nor  broadly,  confidently  and 
cheerfully  proposed. 

We  proceed  to  "give  some  detail  specif  ication  for  this  apparent¬ 
ly  popular  and  promising  prototype  flow. 

.yll!.  The  -prototype  flow. 

The  flow  is  of  axial  symmetry,  and  is  symmetric  fore  and  aft. 

It  can  be  described  in  a  filer idian  plane  using  polar  coordinates  r 
and  (j2  ,  the  origin  of  the  coordinate  system  coinciding  with  the  lump 
center.  The  distance  from  the  axis  of  symmetry  is  then  given  by 
r  sin  (f^ 

The  flow  has  the  radial  component  of  the  velocity  v 


7-i  R  —  'Cos  p  e 

and  a  tangential  velocity  component  at  right  angles  to  the  radius 

7-2  r  =  -  sin  P  (/-£)  £/l 


The  vorticlty  of  the  flow  is  in  consequence 

7-3  g  n  /i  sin  P 

and  thus  its  ratio  to  the  axial  distance  r  isin  £  is 

Je*  (s-2a) 

This  is  constant  over  the  surface  of  all  concentric  spheres 
r  =  const. 

The  axial  velocity  component  is  zero  at  all  points  of  a 


cylinder  having  unit  distance  from  the  axis  of  symmetry.  All  spheres 
about  the  origin  are  surfaces  of  normal  velocity  components  proportional 
to  cos 
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The  flow  shares  this  property  with  the  surface  of  the  Hill  vortex. 

The  momentum  integral  of  the  correlation  flow  converges,  and 
the  total  momehtum  is  zero.  The  specification  for  this  unit  corre¬ 
lation  flow  de-scribes  only  the  type,  but  not  the  intensity.  It 
does  not  indicate  the  specific  kinetic  energy  of  the  turbulence 
flow.  Thai;  depends  on  the  intensity  of  the  correlation  flow,  hamely 
on  the  magnitude  of  its  central  velocity,  at  r  =  0. 

There  are  several  reference  lengths  available  for  the  corre¬ 
lation  flow  of  the  type  (7-1,  (7-2),  the  microscale  of  the 

unit  flow,  Ref.  (t)  is  equal  to  §.  The  square  of  the  radius  of 
the  sphere  dividing  the  points  of  positive  vorticity  from  those  of 
negative  vorticity  is  equal  to  2. 5-  At  the  surface  of  that  sphere, 
the  vorticity  is  zero. 

The  outer  portion  of  the  correlation  flow  merely  describes 
the  statistical  behavior  of  many  adjacent  lumps.  The  inner  portions, 
the  core,  describes  one  average  lump.  The  description  is  particular¬ 
ly  Immediate  near  the  center.  How  far  the  territory  of  one  lump  extends 
depends  on  how  many  lumps  there  are  per  unit  volume.  That  depends 
on  the  (circumstances  of  the  flow  problem.  Broadly  spoken,  the  lumps 
have  a  tendency  to  swell  out  by  viscosity  effects  until  they  come 
so  close  to  each  other  that  the  mutual  viscosity  prevents  further 
swelling  Out,  The  Reynolds  number  of  the  lump  enters  this  picture. 

If  the  Reynolds  number  of  the  lumps  becomes  too  large,  each  lump 
becomes  a  basic  flow  for  smaller  lumps,  This  points  to  an  upper  limit 
of  the  lump  Reynolds  number.  All  this  will  be  discussed  and  computed, 
or  at  least  estimated  in  a  subsequent  paper. 
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IX.  Turbulence  decab. 

Decaying  turbulende 'ife  by  no  means  the  principal  problem  of 
tubulence  reseabcja.  Quasi-steady  turbulence  is  much  more  important. 

The  decaying  turbulence  is  however  always  prominently  featured  in  the 
liteseture  on  the  subject,  because  it  is  the  only  turbulence  problem 
with  respect  to  which  some  meager  and  incomplete  theoretical  results 
have  been  presented.  Experimental  information  in  that  special  field 
is  also  available. 

Only  two  extreme  cases  are  treated  theoretically,  the  initial 
stage  of  turbulence  decay  and  its  final  stage.  The  result  regarding 
the  initial  stage  is  more  acceptable  than  the  other  result. 

The  initial  motion  is  considered  to  be  a  cluster  of  many 
smaller  lumps,  the  cluster  lump  flow  being  the  basic  flow  for  the 
true  lumps .  The  smaller  lumps  take  the  energy  out  of  the  cluster 
and  dissipate  it  soon,  so  that  the  asymptotic  history  of  the  energy 
density  can  be  obtained  from  the  energy  of  the  cluster  or  gross  motion. 
The  clusters  are  thinly  distributed  and  therefore  do  not  interfere  wth 
each  other. 

Under  all  these  circumstances  there  is  an  analogy  between  the 
Velocity  of  an  ordinary  lump  and  the  energy  of  the  cluster  motion. 

Both  are  subject  to  ordinary  diffusion.  The  shear  mixing  m  is  the 
analogue  to  the  kinematic  viscosity  ^  .  It  happens  that  the  final 

result  is  consistent  with  a  constant  mixing  intensity,  in  that  the 
Reynolds  number  of  the  small  lumps  does  not  change.  The  linear  size 
increases  proportional  to  the  square  boot  of  the  time,  and  the  velocity, 
as  we  shall  see,  decreases  inversely  to  the  same. 
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is  not  zero,  the  energy  being  a  positive  quantity,  still  the  rate  of 
decline  of  its  moment  of  inertia  is  small.  It  is  negligible  in  comparison 
with  its  own  magnitude,  because  r  /(t^  )  is  large,  r  denoting  the 
lump  size  and  the  t  the  time  available.  Under  these  circumstances, 
the  moment  of  inertia  of  the  energy  does  not  change  noticeably. 

Considering  it  to  be  dominantly  contributed  by  a  standard  self-pre¬ 
serving  diffusion  distribution  of  finite  moment  of  momentum  of  the 

energy,  the  energy  at  the  cluster  center  would  decline  proportional 
-2.5  -1*5 

to  t  .  But  the  cluster  volume  would  increase  proportional  to  t 

so  that  the  total  energy  per  unit  sja  ce  would  appear  to  be  Inverse 

to  a  linear  function  of  the  time.  This  is  the  result  displayed  in 

the  literature  and  it  seems  to  be  confirmed  experimentally.  There 

are  quite  a  number  of  assumptions  or  limiting  conditions  involved,  but 

they  can  all  be  reasonably  Justified. 

The  opposite  case  is  considered  the  final  stage  of  decay.  The 
current  theory  considers  that  then  viscosity  is  supreme  by  reason 
of  the  small  Reynolds  numbers  prevailing.  That  is  hardly  any  longer 
turbulence,  but  is  of  Interest  as,  the  asymptotic  end  of  turbulence. 

Theory  considers  that  unlimited  space  is  available,  also  unlimited 
time  to  spread  randoihness  all  through  this  Infinite  space.  The  un- 
eveness  of  the  fluid  motion  assumes  .bigger  and  bigger  scale,  without 
limit.  There  regain  always  clusters,,  bigger  and  bigger  in  that  the 
biggest  clusters  are  decaying  at  the1 slowest  rate.  In  accordance 
with  the  reasoning  presented  in  Section  7,  the  energy  decay  will 
therefore  be  finally  Inverse  to  the  2.5th  power  of  the  time. 
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Will  it  really?  What  happens:  finally  depends  entirely  on  the 
boundary  conditions  in  their  widest  sense  including  space  and  time. 

Practically  there  will  always  bdFan  energy  supply  by  way  of  some 
noise  effect.  The  decay  will  then  be  finally  zero.  If  on  the  other 
hand,  the  fluid  is  surrounded  by  hlgid  energy -proof  walls,  the 
decay  will  finally  be  exponential.  Thi  decay  exponent  will  increase 
from  Its  Initial  value  -1  because  the  turbulence  becomes  more  crowded. 
It  is  not  unusual  nor  unreasonable  to.  expect  that  -2.5  will  be  reached 
and  even  surpassed.  The  question  is  whether  the  transient  expondent 
-2.5  is  less  transient  than  the  adjacent  exponents,  so  that  the  decay 
curve  is  in  particularly  intimate  contact  with  a  parabola  correspond¬ 
ing  to  that  exponent. 

If  there  are  good  reasons  to  expect  that,  then  the  finding 
of  the  present  paper  would  make  the  exponent  -2.5  even  more 
acceptable  in  that  the  requirement  for  a  small  Reynolds  number  would 
be  less  severe. 

The  necessary  condition  for  a  good  degree  of  preservation  of 
the  exponent  -2.5  is  that  there  are  closely  packed  lump  clusters 
and  that  there  Is  furthermore  enough  space  and  a  supply  of  even  wider 
spread  velocity  configurations  or  deviations  to  become  the  feeding 
cluster  flows  for  the  less  spread-out  ones. 

The  semi -permanency  of  the  exponent  -2.5  can  not  ba  decisively 
indicated  by  experimental  evidence.  The  precision  of  the  tests  is  not 

t  1 

large  enough.  There  is  however  some  experimental  evidence  to  show 
that  the  exponent  associated  with  the  weakest  energy  intensity 
measureable  is  not  inconsistent  with  an  exponent  of  decay  tfl  =  -2.5- 
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X.  The  ratio  of  the  shear  Intensities. 

Any  statistical  approach  to  the  turbulence  problem  must  have 
as  its  immediate  aim  the  theoretical  determination  of  the  intensities 
of  turbulence  mixing.  This  .mixing  ,  J.6  not  the  primary  subject  of  the 
present  paper.  The  discus&lon  ' of  the  relative  weight  of  the  viscosity 
effects  and  of  the  inertia  effects  must  precede  such  determination. 

The  determination  of  the  intensity  of  mixing  is  a  .subject  which 
requires  its  own  study  and  its  own  paper.  However,  to  offer  in  this 
paper  at  least  some  advance  in  the  direction  of  that  determination, 

'we  will  conclude  it  with  the  discussion  of  the  ratio  of  the  shear 
mixing  to  the  mixing  for  energy  diffusion.  That  these  two  are 
different  from  each  other,  and  why,  will  appear  clearly  in  the  follow¬ 
ing  deduction.  It  will  also  appear  that  the  general  form  used  by 
us  in  Hef .  (m)  is  consistent  with  the  one  immediately  to  be  derived 
for  the  two  mixing  intensities.  The  present  paper  justifies  thus  the 
reasoning  of  this  author  in  hiS  previous  paper  on  the  subject. 

The  shear  mixing  expresses  the  transport  of  flow-wise  momentum 
by  the  lumps  crossing  the  shear  or  gliding  surfaces.  If  the  average 
velocity  of  the  turbulence  flow  is  u,  then 'the  flux  of  the  fluid 
crossing  in  one  or  the  other  directions  is  approximately  u/3 .  The 
momentum  it  transports  depends  On  the  penetration  depth.  This  is 
somewhat  analogous  to  the  mean  free  path  of  the  kinetic  gas  theory, 

i 

but  not  closely  analogous.  We  disbursed  this  Already  in  Ref.  (m) 
and  see  now  clearer  for  what  reason  the  penetration  depth  f_  does  not 
dominantly  depend  on  random  .effects  .but.  .on  vXSoqs ity  effects. 
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The  velocity  of  a  lump  decreases  with  time.  The  rate  of 

-1 

decrease  is  at  least  proportional  to  t,  ,  which  is  for  lumps  spaced 
far  apart  from  each  other.  Ijj  a  quasi-steady  shear  flow,  the  lumps 
are  not  spaced  far  apart,  hut  on  the  contrary,  they  must  be  considered 
closely  packed.  In  that  case,  as  discussed  in  Section  7,  the  reference 

-2.5. 

velocities  decrease  proportional  to  t  TTiis  time,  t,  has  to  be 


counted  from  an  initial  time  point,  such  that  at  the  time  of  the  cross¬ 
ing  t  =  tQ  =  const ( 


where  denotes  any 

reference  length,  Such  as  for  instance  the  microscale. 

i  , 

Hence,  the  perpetration  depth  f  results^proportional  to 

to;  7  •/Tdt  =  d,  A  J0J 


V 

_  M  o 

S* 


Ijh 

v.s-y 


Hence  the  shear  mixing  is  proportional  to 

Mo  =  V 


!0'X 


This  is  exactly  the  form  proposed  in  Hef .  ( m) . 

For  the  diffusion  of  the  energy  we  obtain  a  factor  different 
from  this.  As  the  lump  penetrates  into  a  layer  of  different  energy 
level,  it  does  not  really  transport  its  energy  that  far,  because 
much  of  it  has  been  dissipated  on  the  way.  Hence  an  applicable 
penetration  depth  will  be  obtained,  not  by  integrating  the  velocity, 
but  by  by  integrating  its  product  with  the  energy  still  carried,  and 
dividing  the  outcome  by  the  initial  energy.  Thus 


x- 


-  7.J- 


/0-3 


fit*  E  ~  U°^°  JitJ&J 


6Cp  Up 

3  -  i,  r  y 


A  V<7 
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Hence  it  is  demonstrated  that  the  diffusion  mixing  is  much  smaller 
than  the  shear  mixing.  The  ratio  is  about  3/13. 

The  smaller  diffusion  mixing  would  leave  the  influence  of 
the  secondary  augumentation  pronounced.  However,  the  present 
derivation  is  based  on  the  assumption  that  the  lumps  move  straight 
ahead  on  the  average.  We  have  to  examine  how  far  that  assumption 
is  justifiable  in  shear  flow. 
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